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Summary 

This  study  will  demonstrate  the  importance  of  application  of  CNC  madiines  in  generation 
of  gear  tooth  surfaces  with  new  topology.  This  topology  decreases  gear  vibration  and  will 
extend  the  gear  capacity  and  service  life.  A  preliminary  investigation  by  a  tooth  contact 
analysis(TCA)  program  has  shown  that  gear  tooth  surfaces  in  line  contact  (for  instance, 
involute  helical  gears  with  parallel  axes,  worm-gear  drives  with  cylindrical  worms  etc'.)  are 
very  sensitive  to  angular  errors  of  misalignment  that  cause  edge  contact  and  an  unfavorable 
shape  of  transmission  errors  and  vibration. 

The  new  topology  of  gear  tooth  surfaces  is  based  on  the  localization  of  bearing  contact, 
and  the  synthesis  of  a  predesigned  parabolic  function  of  transmission  errors  that  is  able  to 
absorb  a  piecewise  linear  function  of  transmission  errors  caused  by  gear  misalignment. 

The  report  will  describe  the  following  topics;  (1)  Description  of  kinematics  of  CNC 
machines  with  6  degrees-of-freedom  that  can  be  applied  for  generation  of  gear  tooth  surfaces 
with  new  topology.  (2)  A  new  method  for  grinding  of  gear  tooth  surfaces  by  a  cone  surface 
or  surface  of  revolution  based  on  application  of  CNC  machines.  This  method  provides  an 
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optimal  approximation  of  the  ground  surface  to  the  given  one.  This  method  is  especially 
beneficial  when  undeveloped  ruled  surfaces  are  to  be  ground.  (3)  Execution  of  motions  of 
the  CNC  machine.  The  solution  to  this  problem  can  be  applied  as  well  for  the  transfer  of 
machine-tool  settings  from  a  conventional  generator  to  the  CNC  machine. 

The  developed  theory  required  the  derivation  of  a  modified  equation  of  meshing  based 
on  application  of  the  concept  of  space  curves,  space  curves  represented  on  surfaces,  geodesic 
curvature,  surface  torsion  etc.  Condensed  information  on  these  topics  of  differential  geometry 
is  provided  a.s  well. 


Introduction 

The  design  and  manufacture  of  gears  with  new  topology  of  gear  tooth  surfaces  are  prob¬ 
lems  of  great  importance  for  helicopter  transmissions.  The  existing  technology  of  gears  is 
restricted  with  the  necessity  to  use  cutting  and  grinding  machines  whose  kinematics  is  based 
on  linear  relations  between  the  motions  of  the  tool  and  the  workpiece. 

The  need  of  low-noise  gears  with  increased  load  capacity  and  service  life  can  be  satisfied 
with  a  new  topology  of  gear  tooth  surfaces  that  is  able  to  provide:  (i)  a  reduced  sensitivity  to 
misalignment  and  avoidance  of  edge  contact,  (ii)  a  parabolic  type  of  function  of  transmission 
errors  to  reduce  the  level  of  possible  vibration,  (iii)  a  localized  bearing  contact  with  controlled 
dimensions  of  the  instantaneous  contact  ellipse. 

The  application  of  CNC  (Computer  Numerically  Controlled)  machines  overcomes  the 
obstacles  presented  for  generation  of  gears  with  a  new  surface  topology  by  using  the  existing 
equipment.  The  CNC  machines  are  able  to  provide  computer  controlled  nonlinear  relations 
between  the  motions  of  the  tool  and  the  gear  being  generated.  Although  such  machines  are 
used  at  present  mainly  for  the  installment  of  machine-tool  settings  with  higher  precision, 
their  prosperous  future  is  in  their  application  for  generation  of  gears  with  new  topology.  The 


2 


CNC  machines  are  a  unique  opportunity  for  researchers  to  modify  the  geometry  of  traditional 
gear  drives  and  benefit  industry  with  gear  drives  with  substantially  improved  parameters. 

The  modification  of  geometry  of  gear  tooth  surfaces  requires  from  researchers  a  new 
approach  for  the  development  of  principles  of  conjugation  of  gear  tooth  surfaces.  Application 
of  conjugate  gear  tooth  surfaces  being  in  instantaneous  line  contact  is  in  the  authors’  opinion 
an  anachronism.  Such  gear  tooth  surfaces  are  very  sensitive  to  misalignment  that  causes  the 
shift  of  the  bearing  contact  to  the  edge  and  transmission  errors  of  such  a  type  that  cause 
a  jerk  at  the  transfer  from  one  cycle  of  meshing  to  the  next.  The  considerations  above  are 
true  for  spur  gears,  involute  helical  gears  with  parallel  axes,  and  worm-gear  drives. 

It  is  necessary  as  well  to  change  the  attitude  to  some  finishing  processes  such  as  honing 
and  shaving  applied  for  helical  and  spur  gears.  It  is  not  reasonable  to  require  that  such 
finishing  processes  would  provide  the  exact  screw  involute  surfaces  knowing  ahead  that  only 

modified  tooth  surfaces  are  to  be  applied. 

The  statements  mentioned  above  are  illustrated  with  the  following  drawings. 

Figures  1  to  4  show  the  influence  of  angular  errors  of  misalignment  A7  and  dAp,-  of  involute 
helical  gears  with  parallel  axes  that  cause  edge  contact  (figures  1  and  3)  and  piecewise  almost 
linear  functions  of  transmission  errors  (figures  2  and  4).  The  drawings  above  are  based  on 
the  investigation  performed  by  Reference  [1].  The  design  parameters  of  the  helical  gears  are 
shown  below: 
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Figures  5,  6  and  7  illustrate  the  impact  of  misalignment  of  worm-gear  drives.  The  shift 
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of  the  bearing  contact  is  shown  in  figures  5  and  6,  and  the  undesirable  shift  of  transmission 
errors  is  shown  in  figure  7  that  will  inevitable  cause  vibration.  The  drawings  are  based  on 
the  research  accomplished  by  Reference  [2]. 

The  design  parameters  of  the  worm-gear  drive  are  shown  below: 


iVi 

=  2 

N, 

=  30 

Axial  module,  m 

=  8(mm) 

7 

=  90° 

Shortest  distance,  E, 

=  176(mm) 

Figure  8  illustrates  why  a  predesigned  parabolic  type  of  transmission  errors  is  beneficial 
for  the  gears  with  the  new  topology  (Ref.  [3]).  This  figure  illustrates  the  interaction  of 
a  parabolic  function  of  transmission  errors(provided  at  the  stage  of  synthesis  of  the  gear 
tooth  surfaces)  with  a  linear  function  of  transmission  errors(caused  by  misalignment).  The 
combination  of  these  functions  is  again  a  parabolic  function,  with  the  same  slope  as  the 
predesigned  one  that  is  translated  with  respect  to  the  initial  parabolic  function.  This  means 
that  the  predesigned  parabolic  function  absorbs  the  linear  function  and  keeps  the  shape  of 
a  parabolic  function. 

There  are  three  cases  of  generation  of  the  workpiece  surface  Ep  by  the  given  tool  surface 
Et  by  CNC  machines: 

(1)  Surfaces  Ej  and  Ep  are  in  continuous  tangency,  however  they  contact  each  other  at 
every  instant  at  a  point  not  a  line. 

(2)  Surfaces  Et  and  Ep  are  in  continuous  tangency  and  they  contact  each  other  at  every 
instant  at  a  line.  Surface  Ep  is  generated  in  this  case  as  the  envelope  to  the  family  of 
surfaces  Ej.  The  family  of  surfaces  is  generated  in  relative  motion  of  Et  to  Ep. 

(3)  An  approximate  method  for  generation  of  a  surface  E^  (ground  or  cut)  with  an  optimal 


4 


approximation  to  the  ideal  surface  Sp. 


An  example  of  case  1  is  the  generation,  for  instance,  of  a  die  designed  for  forging  of  a  gear. 
Generation  of  conventional  spiral  bevel  gears  and  hypoid  gears  by  the  Phoenix  machine 
is  the  example  of  case  2  generation.  Case  3  is  the  basic  idea  for  a  new  method  for  surface 
generation  discussed  in  section  4.  Only  cases  2  and  3  of  surface  generation  are  discussed  in 
this  report. 

The  contents  of  the  report  covers  the  following  topics: 

(i)  Description  of  “Phoenix”  and  “Star”  CNC  machines,  that  are  suitable  for  generation 
of  gear  tooth  surfaces  with  new  topology. 

(ii)  Execution  of  motions  of  CNC  machines. 

(iii)  Generation  of  a  surface  with  optimal  approximation  to  the  ideal  surface. 

(iv)  Concept  of  curvatures  that  are  required  for  computations  for  the  proposed  approach 
for  generation. 


1.  “Phoenix”  and  “Star”  CNC  Machines 

“Phoenix”  CNC  Machine 

The  “Phoenix”  CNC  machine  (figure  9)  is  designed  by  the  Gleason  Works  for  generation 
of  spiral  bevel  and  hypoid  gears.  The  machine  is  provided  with  a  total  of  six  degrees-of- 
freedom.  Three  rotational  motions,  and  three  translational  motions  are  used.  The  transla¬ 
tional  motions  are  performed  in  three  mutually  perpendicular  directions.  Two  of  rotational 
motions  are  provided  as  rotation  of  the  workpiece  and  the  rotation  that  enables  to  change 
the  angle  between  the  axes  of  the  workpiece  and  the  tool.  The  sixth  rotational  motion  is 
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provided  as  rotation  of  the  tool  about  its  axis,  and  generally  is  not  related  with  the  pro¬ 
cess  for  generation.  The  motions  with  other  five  degrees-of-freedom  are  provided  as  related 
motions  in  the  process  for  surface  generation. 

Coordinate  Systems  Applied  for  “Phoenix” 

Coordinate  systems  St  (xf,  yt,  Zt)  a^nd  Sp  (xp,  yp,  Zp)  are  rigidly  connected  to  the  tool 
and  the  workpiece,  respectively  (figure  10).  For  further  discussions  we  will  distinguish  four 
reference  frames  designated  in  figure  9  as  /,  II,  III  and  IV.  The  reference  frame  IV  is  the 
fixed  one  to  the  housing  of  the  machine.  Reference  frames  I,  II  and  III  perform  translations 
in  three  mutually  perpendicular  directions,  respectively.  We  designate  coordinate  systems  Sh 
and  Sm.  that  represent  reference  frames  /  and  III,  respectively(figures  9  and  10).  Coordinate 
axes  of  Sh  and  S^,  are  parallel  to  each  other  and  the  location  of  Sh  with  respect  to  S^  is 
represented  by  and  Coordinate  system  St  performs  rotational  motion 

with  respect  to  Sh  about  the  z/i-axis.  To  describe  the  coordinate  transformation  from  Sm  to 
Sp,  we  use  coordinate  systems  Se  and  5(i(figure  10).  Coordinate  system  Se  performs  rotation 
with  respect  to  Sm  about  the  y^-axis.  Coordinate  axes  of  system  Sj  are  parallel  to  the 
respective  axes  of  Se]  the  location  of  origin  Oj  with  respect  to  Og  is  determined  with  the 
parameter  =  const.  Coordinate  system  Sp  performs  rotational  motion  with  respect  to 
Sd  about  the  x^-axis. 

“Star”  CNC  Machine 

A  version  of  the  “Star”  CNC  machine  that  is  provided  with  6  degrees-of-freedom  is  shown 
in  figure  11.  Coordinate  systems  St  (x(,  yt,  Zt),  Sp  {xp,  yp,  Zp)  and  Sj  (x/,  y/,  zy)  are  rigidly 
connected  to  the  tool,  workpiece  and  frame,  respectively.  Coordinate  system  Sd  is  parallel 
to  system  5/  and  the  location  of  Sd  with  respect  to  Sj  is  represented  in  5/  by  (xy^‘'\  0, 
0).  Coordinate  system  Sg  performs  rotational  motion  with  respect  to  Sd  about  the  yd-^xis. 
Coordinate  system  Sh  is  parallel  to  Sg  and  the  location  of  Sh  with  respect  to  Sg  is  represented 
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in  Se  by  (0,  Coordinate  system  St  performs  rotational  motion  with  respect  to 

Sh  about  the  x/^-axis.  Coordinate  system  Sp  performs  rotational  motion  with  respect  to 
the  fixed  coordinate  system  5/  about  the  x/-axis.  Altogether  there  are  three  translational 
motions  along  axes  Xf,  t/e  and  Zg  and  three  rotational  motions  about  axes  x^,  yd,  and  Xh,- 

2.  Basic  Principle  of  Execution  of  Motions  on  CNC  Machine 

Consider  that  the  location  and  orientation  of  the  tool  with  respect  to  the  workpiece  are 
given  in  coordinate  systems  that  are  represented  for  a  conventional  generator  or  for  an  ab¬ 
stract  (mathematical)  model  of  the  process  for  generation.  We  will  consider  for  the  following 
derivations  the  example  of  application  of  the  “Phoenix”  machine.  A  similar  approach  can  be 
applied  for  other  types  of  CNC  machines,  for  instance,  for  the  Star  machine.  Our  goal  is 
to  develop  the  algorithm  for  the  execution  of  motions  of  the  CNC  machine  using  the  initial 
information  mentioned  above.  Reference  [4]  has  used  for  this  purpose  the  existence  of  a 
common  trihedron  for  the  two  couples  of  coordinate  systems  (^St  ,  Sp  and  (^St  ,  Sp  ) 
that  are  applied  for  the  CNC  machine  and  for  the  generating  process,  respectively.  The 
approach  used  in  this  report  is  as  follows: 

(i)  Consider  that  4x4  matrices  mJ^^  and  3  x  3  matrices  (k  =  C,  G)  have  been 
derived.  The  superscripts  “C”  and  “G”  indicate  the  CNC  machine  and  the  abstract 
generating  process,  respectively. 

(ii)  The  matrix  equality 


4?  =  L<f'  (1) 

will  provide  the  same  orientation  of  S^^^  with  respect  to  S^^^  (k  .=  C,  G)  in  both 
reference  frames. 
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(iii)  The  matrix  equality 


Mif’io  0  0  ir=Mi:'io  o  o  i 


V 


(G) 


pi 


will  provide  the  same  position  vector  (OpOt)p  for  both  reference  frames. 

The  application  of  equations  (1)  and  (2)  for  the  execution  of  motions  of  the  “Phoenix” 
machine  is  considered  for  the  two  following  cases:  (i)  a  hypoid  pinion  is  generated  by  appli¬ 
cation  of  a  conventional  generator,  and  (ii)  a  surface  Sp  with  optimal  approximation  to  the 
ideal  surface  Ep  is  generated. 

Derivation  of  Matrix  Lpf^  and  Position  Vector  (OtOp)^*^) 

Using  a  routine  procedure  for  coordinate  transformations,  we  obtain 

cos  n  COS  (f>  —  sin  ^  cos  (p  sin  (f> 

—  cos /i  sin  <;^i  sin  t/j  sin sin  (?!>  sin  ^  cospsinV’ 

=  +  sin  /X  cos  V’  -|-  cos  /x  cos  ib 

—  cos  n  sin  <j)  cos  sin  /x  sin  (f>  cos  ■0  cos  6  cos  0 

—  sin/xsin0  —  cos/xsin0 

We  note  that  hde  and  L„,/,  are  unit  matrices. 

The  derivation  of  the  position  vector  (OjOp)^^^  in  Sp  is  based  on  the  following  consider¬ 
ations: 


(3) 
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(aiio;)^  >  +  (o,Op)f  >  =  {o„o,)<^> 


(i) 


Thus: 

(05;)^  >  =  =  {^)f^  - 

=  x<°-)(i,),  -  xS.°*'(i„),  -  -  4°‘>(k™)p  (4) 

Here:  =  const.,  and  are  considered  as  algebraic  values. 

(ii)  Vector  can  be  represented  in  coordinate  system  with  the  following  matrix 

equation 

10  01"' 

-i/!?‘'Lp™[0  1  0  f  -  z>?*>L,„[  0  0  (5) 

where  =  LpjLj.L,„  (Lj,  is  a  unitary  matrix). 

Equation  (5)  yields 

x^*^**)  —  cos  <f>  —  z^'*^  sin  <f> 

x^'')  sin  ^  sin  V’  -  cos  ^  -  z^'*)  cos  (f)  sin  V’  (^) 

sin  ^  cos  V’  +  sin  V’  -  z^'*^  cos  (j)  cos  ip 
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3.  Example:  Generation  of  Hypoid  Pinion  by  “Phoenix” 


Generation  of  Pinion  Tooth  Surface  by  Conventional  Generator 

The  pinion  tooth  surface  is  generated  as  the  envelope  to  the  family  of  tool  surfaces  that 
are  cone  surfaces  as  shown  in  figure  12. 

Henceforth,  we  will  consider  the  following  coordinate  systems:  (i)  the  fixed  ones,  So  and 
Sq  that  are  rigidly  connected  to  the  cutting  machine  (figures  13  and  14);  (ii)  the  movable 
coordinate  systems  Sc  and  Sp  that  are  rigidly  connected  to  the  cradle  of  cutting  machine 
and  the  pinion,  respectively;  (iii)  coordinate  system  St  that  is  rigidly  connected  to  the  head 
cutter.  In  the  process  of  generation  the  cradle  with  Sc  performs  rotational  motion  about  the 
2o-axis  with  angular  velocity  and  the  pinion  with  Sp  performs  rotational  motion  about 

the  I, -axis  with  angular  velocity  (figure  14). 

The  tool  (head-cutter)  is  mounted  on  the  cradle  and  performs  rotational  motion  with  the 
cradle.  Coordinate  system  St  is  rigidly  connected  to  the  cradle.  To  describe  the  installment 
of  the  tool  with  respect  to  the  cradle  we  use  coordinate  system  Si>  (figures  12  and  13).  The 
required  orientation  of  the  head-cutter  with  respect  to  the  cradle  is  accomplished  as  follows: 
(i)  coordinate  systems  Sb  and  St  are  rigidly  connected  and  then  they  are  turned  as  one  rigid 
body  about  the  z^-axis  through  the  swivel  angle  j  =  2Tr  —  6  (figure  13);  (ii)  then  the  head- 
cutter  with  coordinate  system  St  is  tilted  about  the  3/i,-axis  under  the  angle  i  (figure  12(b)). 
The  head-cutter  is  rotated  about  its  axis  Zt  but  the  angular  velocity  in  this  motion  is  not 
related  with  the  generation  process  and  depends  only  on  the  desired  velocity  of  cutting. 

The  pinion  setting  parameters  are  Em-  the  machine  offset,  jm-  the  machine-root  angle, 
AB-  the  sliding  base,  and  AA-  the  machine  center  to  back  are  shown  in  figure  14.  The 
head- cutter  settings  parameters  are  Sr-  radial  setting,  6c-  initial  value  of  cradle  angle,  j- 
the  swivel  angle(figure  13),  and  i-  the  tilt  angle(figure  12(b)). 

Pinion  Tool  Surface  Equations 
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The  head-cutter  surface  is  a  cone  and  is  represented  in  St  (figure  12(a))  as 

(tc  +  s  sin  a)  cos  9 
(tc  +  s  sin  a)  sin  9 

rt{s,9)  = 

—s  cos  a 

1 

Here:  (s,  9)  are  the  Gaussian  coordinates,  a  is  the  blade  angle  and  is  the  cutter  point 
radius.  Vector  function  (7)  with  a  positive  and  a  negative  represents  surfaces  of  two  head- 
cutters  that  are  used  to  cut  the  pinion  concave  side  and  convex  side,  respectively. 

The  unit  normal  to  the  head-cutter  surface  is  represented  in  St  by  the  equations 

nt  =  [  —  cos  a  cos  9  —  cos  a  sin  9  —  sin  a  (8) 

The  family  of  tool  surfaces  is  represented  in  Sp  by  the  matrix  equation 

rp(s,  9,  (f)p)  =  Mpg  Mgn  M„o  Moc  Met  Mtt  rt(s,  9)  (9) 

Here:  5n  is  an  auxiliary  fixed  coordinate  system  whose  axes  parallel  to  So  axes. 

cosz  0  sinz  0 

0  10  0 

—  sin  z  0  cos  z  0 

0  0  0  1 
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—  sinj  —cosj  0  Sr 

cosj  —sinj  0  0 

= 

0  0  10 
0  0  0  1 

cos  q  sin  g  0  0 

—  sin  q  cos  5  0  0 

0  0  10 

0  0  0  1 

10  0  0 

0  1  0 

0  0  1  -AS 
0  0  0  1 

cos7„  0  sin  7^  -AA 

0  10  0 

-sin7„i  0  cos7„  0 

0  0  0  1 

10  0  0 

0  cos  <j>r,  —  sin  (j)-p  0 

IVI  = 

iVipj 

0  sin  (f>j,  cos  (f>p  0 

0  0  0  1 

6  =  2t  —  j]  q  =  6^  -{■  mcp(f)p  where  6^  is  the  initial  cradle  angle  and  m^p  = 

Equation  of  Meshing 

This  equation  is  described  in  Reference  [5]  as: 

n(p)  .  =  f{s,  e,  4>p)  =  0  (10) 
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where  n^'^and  'are  the  unit  normal  and  the  normal  to  the  tool  surface,  and  is  the 
relative  velocity  between  the  tool  surface  and  the  workpiece. 

Equation  (10)  is  invariant  with  respect  to  the  coordinate  system  where  the  vectors  of  the 
scalar  product  are  represented.  These  vectors  in  our  derivations  have  been  represented  in  So 
as  follows, 


iio  —  Loc  Eci  Ljt  rif 


=  [(  -  wl"))  X  r,]  +  {OoA  X 


Here: 


To  =  Moc  Mci  Mbt  Tt 

0:4  =  [0  -Em  AJ5]^ 

=  -[  cos 7  0  sin7  ;  (|  =  1) 

a;W  =  -[00  mop  f 
Pinion  Tooth  Surface 

Equations  (9)  and  (10)  represent  the  pinion  tooth  surface  in  three-parametric  form  with 
parameters  s,  9  and  ^p.  However,  since  equation  (10)  is  linear  with  respect  to  s  we  can 
eliminate  s  and  represent  the  pinion  tooth  surface  in  two-parametric  form  as 

rp{e,<f>p,dk)  (11) 
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Here:  dk  {k  =  1, . . .  ,8)  designate  the  installment  parameters:  Em,  7m,  AB,  AA,  Sr,  9c, 
j  and  i. 

The  normal  to  the  pinion  tooth  surface  is  represented  as 

np{6,(f>p,dk)  (12) 

where  dk  {k  =  1,2, 3, 4)  designate  the  installment  parameters  7„,  6c,  j  and  i. 

Derivation  of  and  {OtOpYJ^^ 

Our  next  goal  is  to  derive  the  algorithm  for  execution  of  motions  on  “Phoenix  ” ,  knowing 
the  basic  machine-tool  settings  on  the  conventional  generator. 

The  coordinate  systems  applied  for  the  CNC  machine  are  represented  in  figure  10.  The 
performed  coordinate  transformation  yields: 

=  Mq)]  (fc  =  1, 2, 3;  /  =  1, 2, 3)  (13) 

Here: 

an  =  cos  i  cos^m  sin(g  -  j)  -  sin  i  sin  7m  , 
an  =  -  cos(9  -  ;■)  cosjm  , 
an  =  sin  i  cos  7m  s\n{q  —  j)  +  cos  i  sin  7m  , 

021  =  cos  i  sin  7m  sin  4>p  sin(5  —  j)  +  cos  i  cos{q  —  j)  cos  (f>p 
-f  sin  i  cos  7m  sin  (t>p  , 

022  =  -  cos(g  -  j)  sin  7m  sin<i!>p  +  sin(5  -  j)  cos  4>p  , 

023  =  sin  i  sin  ~fm  sin  sin(g  —  j)  +  sin  i  cos[q  —  j)  cos  (f>p 

—  cos  i  cos  7m  sin  (f)p  , 

031  =  —  COS  i  sin  7m  cos  <f)p  sm{q  —  j)  +  cos  i  cos{q  —  j)  sin  <i>p 

—  sin  i  cos  7m  cos  (f)p  , 

032  =  sin  7m  cos  (f)p  cos{q  -  j)  +  sm{q  -  j)  sin  <t>p  , 

033  =  —  sini  sin  7m  sm{q  —  j)  cos  6p  +  sini  cos(g  —  j)  sin  </>p 
-f  cos  i  cos  7m  cos  <f)p 
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The  variable  parameters  q  and  (^p  are  related  and  therefore  the  coefficients  a),i  {k  = 
1,2,3;/  =  1,2,3)  are  functions  of  q. 

The  position  vector  (OtOp)p^^  is  represented  as  follows: 


(OtOp)f )  =  0  0  0  1 

Sr  cos  q  cos  —  AjB  sin  7^  —  A  A 
— 5'ii(sin  q  cos  ^p  —  cos  q  sin  'y^  sin  (j>p) 
+Em  cos  (f>p  +  AjB  cos  7m  sin  (f>p 
— 5fl(sin  q  sin  (f>p  +  cos  q  sin  7^  cos  <j>p) 
+E^  sin  ^p  -  AB  cos  7^  cos  (f>p 
1 

014(9) 

_  024(9) 

034(9) 

1 


(15) 


Execution  of  Motions  of  CNC 

Matrix  equality  (equation  (1))  provides  nine  dependent  equations  for  determination  of 
functions  <l){q),  il^{q),  and  fi{q).  We  can  determine  these  functions  using  the  following  proce¬ 
dure: 

Step  1:  Determination  of  (f>. 


sin  <f>  =  013(5) 


(16) 


This  equation  provides  two  solutions  (f>]  the  smaller  value  of  (j)  can  be  chosen. 
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step  2:  Determination  of  V’- 


cos  ^  sin  V’  =  ^23(9)  5  cos  ^  cos  ^  =  033(5)  (17) 

These  equations  provide  a  unique  solution  for  V’,  considering  ({)  as  given. 

Step  3:  Determination  of  //. 

cos /i  cos  ^  =  011(5)  ,  —  sin// cos  (?i>  =  012(5)  (18) 

These  equations  provide  a  unique  solution  for  //,  considering  4>  as  given. 

For  the  generation  a  face-milled  hypoid  pinion,  a  tool  with  a  cone  surface  is  applied.  The 
tool  surface  is  a  surface  of  revolution  and  the  rotation  of  the  tool  about  its  axis  is  not  related 
with  4>.  Functions  (17)  must  be  applied  and  executed  only  for  the  generation  of  face-hobbed 
hypoid  pinion,  that  is  cut  by  a  blade. 

Vector  equality 


(0,0,)fl  =  (0,0,)f'  (19) 

permits  the  determination  of  functions  a:|^'‘^(5),  Equations  (6),  (15) 

and  (19)  considered  simultaneously,  represent  a  system  of  three  linear  equations  in  the  un¬ 
knowns:  The  solution  to  these  equations  enables  to  determine  the  trans¬ 

lational  motions  on  the  CNC  machine. 

4.  Generation  of  a  Surface  with  Optimal  Approximation  To  the  Ideal  Surface 

Introduction 

This  section  is  based  on  the  research  accomplished  by  Reference  [2]  that  was  directed  at 
generation  of  a  surface  (S^)  that  must  be  in  optimal  approximation  to  the  theoretical  (ideal) 
surface  Sp. 
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.  The  method  for  generation  of  Eg  is  based  on  following  ideas: 

(1)  A  mean  line  Lm  on  the  ideal  surface  Sp  is  chosen  as  shown  in  figure  15. 

(2)  The  tool  surface  St  is  a  properly  designed  surface  of  revolution  (in  particular  cases  St 
is  a  circular  cone  as  shown  in  figure  15)  that  moves  along  Lm-  Surfaces  St  and  Sp  are 
in  continuous  tangency  along  Lm'i  M  is  the  current  point  of  tangency  (figure  15).  The 
orientation  of  S*  with  respect  to  Sp  (determined  with  angle  p)  is  continuously  varying. 
Angle  13  at  current  point  M  of  tangency  is  formed  by  the  tangents  t/  and  tj  to  Lm  and 
the  tool  generatrix,  respectively  (figure  15).  Tangents  tf  and  form  plane  11  that  is 
tangent  to  St  and  Sp  at  point  M. 

(3)  The  tool  surface  St  in  its  motion  with  respect  to  Sp  swept  out  a  region  of  space  as  a 
family  of  surfaces  St.  The  envelope  to  the  family  of  St  is  surface  Sg,  the  ground  or 
cut  surface,  that  is  in  tangency  with  the  theoretical  surface  Sp  at  any  point  M  of  Lm 
and  must  be  in  optimal  approximation  to  Sp  in  any  direction  that  differs  from  Lm- 

(4)  The  optimal  approximation  of  Sg  to  Sp  is  obtained  by  variation  of  angle  jS  (figure  15). 

(5)  The  continuous  tangency  of  tool  surface  Sf  with  Sp  and  properly  varied  orientation 
of  St  can  be  obtained  by  the  execution  of  required  motions  of  the  tool  by  a  computer 
controlled  multi-degree-of-freedom  machine.  One  of  these  degrees  of  freedom,  rotation 
of  the  tool  about  its  axis,  provides  the  desired  velocity  of  grinding  (cutting)  and  is  not 
related  with  the  process  for  generation  of  Sg. 

The  contents  of  this  section  cover  the  following  topics: 

(1)  Determination  of  equation  of  meshing  between  the  tool  surface  S*  and  the  generated 
surface  Sg.  The  equation  of  meshing  provides  the  necessary  condition  of  envelope 
existence  to  the  family  of  surfaces. 


17 


(2)  Determination  of  generated  surface  as  the  envelope  to  the  family  of  surfaces  Et 

swept  out  by  the  tool.  Surface  coincides  with  the  theoretical  (ideal)  surface  Sp 

along  the  mean  line  and  deviates  from  Ep  out  of  Lm. 

(3)  Determination  of  deviations  of  Eg  from  Ep  (in  regions  that  differ  from  Lm)  and  mini¬ 
mizations  of  Eg  deviations  for  optimal  approximation  of  Eg  to  Ep. 

(4)  Determination  of  curvatures  of  Eg  that  are  required  when  the  simulation  of  meshing 
and  contact  of  two  mating  surfaces  are  considered. 

(5)  Execution  of  required  motions  of  Et  with  respect  to  Ep  by  appHcation  of  a  multi-degree- 
freedom,  computer  numerically  controlled  machine. 

An  effective  approach  for  the  derivation  of  the  necessary  condition  of  the  envelope  Eg 
existence  is  discussed.  This  method  is  based  on  the  idea  of  motion  of  the  Darboux-Frenet 
trihedron  along  Lm,  the  chosen  mean  line  of  Ep. 

An  additional  effective  approach  for  determination  of  curvatures  of  generated  surface  Eg 
is  discussed  as  well.  This  approach  is  based  on  the  fact  that  the  normal  curvatures  and 
surface  torsions  (geodesic  torsions)  of  Eg  are:  (i)  equal  to  the  normal  curvatures  and  surface 
torsions  of  Ep  along  Lm,  and  (ii)  equal  to  the  normal  curvatures  and  surface  torsions  of  tool 
surface  Ej  along  the  characteristic  Lg  (the  instantaneous  line  of  tangency  of  Ej  and  Eg). 

Mean  Line  on  Ideal  Surface  Ep 

The  ideal  surface  Ep  is  considered  as  a  regular  one  and  is  represented  as 

r,(u„9J€C^  (20) 

where  (up,  9p)  are  the  Gaussian  coordinates  of  Ep. 

The  unit  normal  to  Ep  is  represented  as 
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dvp  dvp 

^p  =  a  ^  ^ 
dup  oOp 


(21) 


The  determination 


N. 


np  = 


INpl  ’ 

of  mean  line  on  L-m 


is  based  on  the  following  procedure: 


(i)  Initially,  we  determine  numerically  n  points  on  surface  Sp  that  will  belong  approxi¬ 
mately  to  the  desired  mean  line  Lm.- 

(ii)  Then,  we  can  derive  a  polynomial  function 


i=i 


{i  =  l,...,n) 


(22) 


that  will  relate  surface  parameters  (up,^p)  for  the  n  points  of  the  mean  line  on  Sp. 


The  mean  line  L,„,  tangent  Tp  and  unit  tangent  tp  to  the  mean  line  are  represented  as 
follows 


^p(^p(^p)’ ^p)  ’ 


dvp  dup 
dup  dOp 


tp 


(23) 


The  constraint  for  tp  is  that  it  must  be  of  the  same  sign  and  differ  from  zero  at  the  same 
intervals  of  interpolation. 

Tool  Surface 

The  tool  surface  St  is  represented  in  coordinate  system  St  rigidly  connected  to  the  tool 
by  the  following  equations 


Xf  =  Xt(ut)  cos  9t  t  Vt  —  sin  9t  ,  Zt  —  zt{ut) 


(24) 
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The  axial  section  of  Ej  obtained  by  taking  =  0  represents  a  circular  arc,  or  a  straight 
line  in  the  case  when  E*  is  a  circular  cone.  Surface  as  shown  in  equations  (24)  of  the  tool  is 
formed  by  rotation  of  the  axial  section  of  E(  about  the  z^-axis. 

The  surface  unit  normal  is  determined  as 

Ni  dvt  dvt 

|Nf|  odf  oui 

Equation  of  Meshing  Between  Et  and  E^ 

Equation  of  meshing  represents  the  necessary  condition  of  existence  of  envelope  E^  to 
the  family  of  surfaces  Et  that  is  swept  out  by  the  tool  surface  E^. 

The  equation  of  meshing  can  be  derived  by  using  the  equation 

=  0  (26) 

Here:  i  indicates  the  coordinate  system  E,-  where  the  vectors  of  the  scalar  product  are 
represented;  is  the  normal  to  surface  E*;  is  the  relative  velocity  in  the  motion  of 
E(  with  respect  to  E^. 

Henceforth,  we  will  consider  two  basic  coordinate  systems,  St  and  5p,  that  are  rigidly 
connected  to  the  tool  surface  Et  and  the  ideal  surface  Ep.  In  addition  to  Et,  we  will  consider 
two  trihedrons:  5i,(tf,, dj, nj,)  and  5/(ty^, dy, n^).  Trihedron  Sj,  is  rigidly  connected  to  Et  and 
coordinate  system  St  (figure  16).  Here:  Oh  is  the  point  of  the  chosen  generatrix  of  Et  where 
the  trihedron  is  located;  tj,  is  the  tangent  to  the  generatrix  at  Ob',  n;,  is  the  surface  unit 
normal  of  Et  at  Ob',  dj,  =  nj,  x  t{,;  vectors  t;,  and  d{,  form  the  tangent  plane  to  Et  at  Ob- 
Trihedron  Sf  moves  along  the  mean  line  (figure  17);  ty  is  the  tangent  to  the  mean  line 
at  current  point  M  (figure  17);  ny  is  the  surface  unit  normal  to  Ep  at  point  M;dy  =  nyxty; 
vectors  ty  and  dy  form  the  tangent  plane  to  Ep  at  point  M. 

The  tool  with  surface  Et  and  trihedron  St  moves  along  mean  line  Lm  of  Ep  and  Ob 


(25) 
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coincides  with,  current  point  M  of  mean  line  L-m-  Surfaces  St  and  Sp  are  in  tangency  at  any 
current  point  M  of  mean  line  L^.  The  orientation  of  Sh  with  respect  to  St  is  determined 
with  angle  ^  that  is  varied  in  the  process  for  generation. 

We  start  the  derivations  with  the  case  when  St  is  a  circular  cone(figure  18).  The  angular 
velocity  w/  of  rotation  of  5/  with  respect  to  is  represented  as 

ds 

=  (27) 

Here:  i  is  the  surface  torsion  (geodesic  torsion),  and  kg  are  the  normal  and  geodesic  cur¬ 
vatures  of  surface  Sp  at  current  point  M  of  mean  line  L^.,  ds  is  the  infinitesimal  displacement 
along  Lm- 

The  angular  velocity  Of  of  trihedron  Sb  is  represented  in  Sf  as 

dd  dB  ds 

+ 'S'"' =  I ‘ 


The  orientation  of  cone  Sf  is  determined  by  function  Bi&p)  and 

^  _  dB.  1 

ds  dOp  ds  ddp  |Tp| 


(29) 


where  Tp  is  the  tangent  to  the  mean  line  at  current  point  M. 

The  transformations  of  vector  components  in  transition  from  St  to  Sb  and  Sb  to  Sf  are 
represented  by  3  x  3  matrix  operators  Ltt  and  L/j,.  Here: 

-  sin  B  0 
cos  B  0 
0  1 


cos  B 
L/t  =  sin^ 
0 
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sin  7t  cos  6t  sin  7t  sin  6t  cos  7* 

L(,f  =  sin^t  —  cos^t  0 

cos  7t  cos  6t  cos  7t  sin  6t  —  sin  7t 

The  cone  surface  Ht  is  represented  in  St  as  follows  (fig.  18) 

■iT 

ft  =  Ut  sin7tCos^t  sin  7*  sin  cos  74  (32) 

where  {ut,  0<)  are  the  surface  parameters,  74  is  the  cone  apex  angle. 

The  unit  normal  to  the  cone  surface  is 

-iT 

n4  =  cos  74  cos  Of  cos  74  sin  Of  —  sin  74  (33) 

The  sought-for  equation  of  meshing,  necessary  condition  of  existence  of  envelope  is 
represented  in  the  form; 


nf  ■  =  0 


where 


=  L/4n4  (35) 

The  derivation  of  expression  is  simplified  while  taking  into  account  the  following 
considerations: 

(a)  The  relative  velocity  vector  can  be  represented  as 


V/  -  Tf  -t  Zf 
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Here,  ^  is  the  skew- symmetric  matrix  represented  as 


0  — W3  1-02 

Ct>3  0  — U}\ 


— (jJ2  0 


Vector  fif\s  represented  by 


(37) 


f2 f  =  Wit f  LJ2<if  LOzUf  — 


t  -K 


diYd^ 

ds  \  it 


(38) 


(b)  Consider  that  point  N  on  surface  St(fig.  15)  is  the  point  of  the  characteristic  (the 
line  of  tangency  of  St  and  the  generated  surface  S^).  Certainly,  the  equation  of  meshing 
must  be  satisfied  for  point  N. 

The  position  vector  OfN  can  be  represented  as(figs.  15  and  18) 


OfN  =  OtN  -  OtOf  (39) 

Here,  0^  is  the  position  vector  of  point  N  that  is  drawn  from  the  origin  Ot  of  St  to  N-, 


vector  OtN  is  represented  in  St  as 


OtN  =  UtBt  =  Ut(sin7t  cos^f  if  -1- sin  7f  sin  ^t  jf  -|-  cos7f  kf) 


(40) 


where 


Gf  = 


(41) 


is  the  unit  vector  of  cone  generatrix  OtN. 

Vector  OtOf  (figure  18)  is  represented  in  Sb  as 
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OtOf  —  It  4 


(42) 


where  /<  =  \OtOf 


Vector  OfN  is  represented  in  5/  using  the  matrix  equation 


—  UtTift^t  ~  h^fbh 

(c)  We  represent  now  the  equation  of  meshing  as 


(43) 


n<;> .  v<“>  =  nW  .  -  (,Ly,ij)]  +  (nf  ■  t,) 


ds 

dt 


(44) 


(d)  The  further  simplification  of  equation  of  meshing  is  based  on  the  following  rule  for 
operations  with  skew-symmetric  matrices  [5]: 


(45) 


Here:  and  designate  skew- symmetric  matrices,  A^  is  the  transpose  matrix  for  A. 

Considering  that  elements  of  B^®^  are  represented  in  terms  of  components  of  vector 


b  =  [  6i  b2  63  j  (46) 

we  obtain  that  the  elements  of  skew- symmetric  matrix  C^®^  are  represented  in  terms  of 
components  of  vector  c,  where 


r 

T 

Cl  C2  C3 

=  aJ 

61  62  63 

iT 


(47) 


ds 


Using  the  above  considerations  and  eliminating  — ,  the  final  expression  of  equation  of 

dt 

meshing  can  be  represented  as 
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Here: 


n^P  ■  =  f{ut,  6t,  6p)  =  utnf  +  nf  —  0 


(48) 


A<-.^  =  LJ.«<*>L 


ft 


B(*) 


ds 

dt 


‘fb 


(49) 


0 

—03 

02 

as 

0 

-Oi 

-02 

Oi 

0 

(50) 


(51) 


0 

-h 

62 

— t  cos  +  fcn  sin  ^ 

b(")  = 

^3 

-h 

0 

6i 

-w 

0 

62 

63 

— 

t  sin  13  +  kn  cos  /3 

.  +  . 

The  family  of  characteristics  Lg  and  the  instantaneous  lines  of  tangency  of  Sj  and  are 
represented  in  St  by  the  equations 


Cl 

02 

as 


d/3. 

t  cos  ^  sin  7t  —  kn  sin  ^  sin  jt  +  [kg  +  —)  cos  jt 


tsin^  +  kn  cos  ^ 


d/3 

t  cos  ^  cos  7t  —  kn  sin  ^  cos  Jt  —  (kg  +  — )  sin  7t 


Ft  =  rt(ut,0t)  ,  fiuu^t.Op)  =  0  (53) 

where  9p  is  the  parameter  of  the  family  of  Lg.  Taking  6p  =  0^  (i  =  l,2,...,n.),  we  obtain 
the  current  characteristics  on  surface  St. 
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It  is  easy  to  verify  that  the  equation  of  meshing  between  and  T,g  is  satisfied  the  current 
point  M  of  mean  line  Lm  on  the  ideal  surface  Ep.  This  means  that  the  characteristic  Lg 
intersects  Lm  at  point  M,  for  which  we  can  take  =  0  since  is  a  surface  of  revolution.  In 


the  case  when  Sj  is  a  circular  cone  (figure  18),  we  can  take  for  point  M  that  Ut  =  =  /<. 

The  approach  discussed  above  for  the  derivation  of  the  equation  of  meshing  can  be  easily 
extended  for  application  in  the  more  general  case  when  the  tool  surface  is  a  general  surface 
of  revolution. 


Determination  of  Generated  Surface 

The  ground  surface  is  generated  as  the  envelope  to  the  family  of  tool  surface  S(; 
surface  E^  is  represented  in  5p  by  the  following  equations 


=  o  (54) 

Here:  f{ut, =  0  is  the  equation  of  meshing;  is  the  equation  of  tool  surface 

Sf  represented  in  5/;  r^^^(iip(0p), 0p)  is  the  vector  function  that  represents  in  Sp  the  mean 
line  Lm',  the  3x3  matrix  operator  Lp/  which  transforms  vectors  in  transition  from  Sf  to  Sp 
is  represented  as 


Lp/  — 


■px 

d-px 

^px 

t 

•py 

dpy 

Upy 

t 

■pz 

dpz 

ripz 

where 


is  the  unit  tangent  to  the  mean  line  Lm', 


(55) 


(56) 
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np 


=  ± 


dvp  dtp 
dup  ^  dOp 


dTp 

dup  d6p 


(57) 


dp  =  Hp  X  tp  (58) 

The  sign  chosen  in  equation  (57)  must  provide  the  direction  of  np  toward  to  the  surface 
of  the  workpiece  under  consideration. 

Equations  (54)  represent  in  Sp  the  generated  surface  Sp  in  three-parametric  form  but 
with  related  parameters.  Parameter  ut  is  linear  in  equation  of  meshing  when  Et  is  a  cone, 
therefore  this  parameter  can  be  eliminated  and  the  generated  surface  Sp  can  be  represented 

in  Sp  as 


=  Tff  =  rs(^P.^t)  (59) 

Remember  that  surfaces  Sp  and  Sp  have  a  common  line  Lm  where  they  are  in  tangency. 
Surface  Sp  is  in  tangency  with  Sj  along  the  instantaneous  line  Lg  that  passes  through  current 
point  M  of  line  L^.  The  tangents  to  Lg  and  Lm  lie  in  the  plane  that  passes  through  M  and 
is  tangent  to  three  surfaces  (Sp,Sp  and  St)  simultaneously. 

Optimal  Approximation  of  Generated  Surface  Sp  to  Ideal  Surface  Sp 

The  procedure  of  optimal  approximation  of  Sp  to  Sp  is  divided  into  the  following  stages: 
(i)  design  of  grid  on  Sp,  the  net  of  points,  where  the  deviation  of  Sp  from  Sp  will  be 
determined;  (ii)  determination  of  initial  function  for  the  first  iteration;  angle  /? 

determines  the  orientation  of  the  tool  surface  Sj  with  respect  to  Sp  (figures  15  and  17),  (iii) 
determination  of  deviations  of  Sg  from  Sp  with  the  initial  function  /^'^H^p);  (iv)  optimal 
minimization  of  deviations. 

Grid  on  Surface  Sp.  Figure  19(a)  shows  the  grid  on  surface  Sp,  the  net  of  (n,m)  points, 
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where  the  deviations  of  Sj  from  Sp  are  considered.  The  position  vector  is  OpQij  = 
(figure  19(b)).  The  computation  is  based  on  the  following  procedure: 

(i)  The  desired  components  Lij  and  Rij  of  the  position  vector  are  considered  as 
known. 

(ii)  Taking  into  account  that 


Li,,  =  ,  Rli  =  ixW'K,  e,)f  +  e,)Y 


(60) 


we  will  obtain  the  surface  Ep  parameters  for  each  grid  point. 


Determination  of  Initial  Function  (^p)-  The  determination  of  is  based  on 

the  following  idea:  the  instantaneous  direction  of  tj  (the  tool  generatrix)  with  respect  to 
tangent  t/  to  the  mean  line  Lm  (figure  17)  must  provide  the  minimal  value  Here: 

is  the  relative  normal  curvature  determined  eis 


-  jbip) 


(61) 


where  and  are  the  normal  curvatures  of  surfaces  Sj  and  Sp  along  tj,.  In  the  case  of 
nondevelopable  ruled  surface  Ep,  vector  t;,  can  be  directed  along  the  asymptote  of  Sp. 

The  requirement  that  is  minimal,  enables  to  determine  function  ^^^\6p)  numerically. 
Since  we  need  for  further  computations  the  derivative  function  ^^^\6p)  is  represented 

dvp 

as  a  polynomial  function  that  must  satisfy  the  numerical  data  obtained  for  the  chosen  points 
of  mean  line  L^. 

Determination  of  Deviations  of  E^  from  Ep.  We  are  able  at  this  stage  of  investigation  to 
determine  the  equation  of  meshing  between  surfaces  Sf  and  Sp,  and  surface  Sp  as  discussed 
above.  The  computation  of  deviations  of  Sp  from  Sp  at  the  grid  points  is  based  on  the 
following  considerations: 


28 


(i)  Surfaces  Ep  and  are  represented  in  the  same  coordinate  system  (5p)  by  the  following 
vector  functions: 


rp{up,6p)  ,  rg{9g,et)  (62) 

(ii)  The  position  vector  and  surface  coordinates  are  known  for  each  point 

Q(».i)  of  the  grid  on  surface  Ep. 

(iii)  Point  on  surface  Ep  corresponds  to  point  on  surface  Ep.  The  surface  Ep 

parameters  can  be  determined  by  using  the  following  two  equations 


(63) 


(iv)  Due  to  deviations  of  S,  from  Ep,  we  have  that  ^  i6.i).  The  deviation  of  E,  from 
Ep  at  the  grid  point  is  determined  by  the  equation 


«i,,  =  n<«)-(rW)-rW>)  (64) 

where  nj,'"’*  is  the  unit  normal  to  surface  Ep  at  the  grid  point 

The  deviation  Sij  can  be  positive  or  negative.  We  designate  as  positive  such  a  deviation 
when  6i^j  >  0  considering  that  is  directed  into  the  “body”  of  surface  Ep.  Positive  devi¬ 
ations  of  Ep  with  respect  to  Ep  provide  that  Ep  is  inside  of  Ep  and  surface  Ep  is  “crowned”. 

It  is  not  excluded  that  initially  the  inequality  >  0  is  not  observed  yet  for  all  points 
of  the  grid.  Positive  deviations  6ij  can  be  provided  choosing  the  following  options: 

(1)  choosing  a  surface  of  revolution  with  a  circular  arc  in  the  axial  section  instead  of  a  circular 
cone;  a  proper  radius  of  the  circular  arc  must  be  determined. 

(2)  changing  parameter  =  \OtOb\  (figures  17  and  18);  this  means  that  the  grinding  cone 
will  be  displaced  along  with  respect  to  the  mean  line  Lm- 
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(3)  varying  the  initially  chosen  function 

Minimization  of  Deviations  Consider  that  deviations  Sij  {i  =  j  = 

of  Ep  with  respect  to  Ep  have  been  determined  at  the  {n,m)  grid  points.  The  minimization 
of  deviations  can  be  obtained  by  corrections  of  previously  obtained  function  The 

correction  of  angle  is  equivalent  to  the  correction  of  the  angle  that  is  formed  by  the  principal 
directions  on  surfaces  Ej  and  E^.  The  correction  of  angle  ^  can  be  achieved  by  turning  of 
the  tool  about  the  corrunon  normal  to  surfaces  Sf  and  Ep  at  their  instantaneous  point  of 
tangency  Mk- 

The  minimization  of  deviations  is  based  on  the  following  procedure; 

Step  1:  Consider  the  characteristic  Lgk,  the  line  of  contact  between  surfaces  E(  and  Ep, 
that  passes  through  current  point  Mk  of  mean  line  on  surface  Ep  (figure  20).  Determine 
the  deviations  Sk  between  Ej  and  Ep  along  line  Lgk  and  find  out  the  maximum  deviations 
designated  as  Points  of  Lgk  where  the  deviations  are  a  maximum  are 

designated  as  and  These  points  are  determined  in  regions  I  and  II  of  surface  Ep 
with  line  Lm  as  the  border.  The  simultaneous  consideration  of  the  maximum  deviations  in 
both  regions  permits  the  minimization  of  the  deviations  for  the  whole  surface  Ep. 

Note:  The  deviations  of  E(  from  Ep  along  Lgk  are  simultaneously  the  deviations  of  Ep 
from  Ep  along  Lgk  since  Lgk  is  the  line  of  tangency  of  Ef  and  Ep. 

Step  2:  The  minimization  of  deviations  is  accomplished  by  correction  of  angle  that  is 
determined  at  point  Mk  (figure  20).  The  minimization  of  deviations  is  performed  locally, 
for  a  piece  k  of  surface  Ep  with  the  characteristic  Lgk-  The  process  of  minimization  is  a 
computerized  iterative  process  based  on  the  following  considerations: 

(i)  The  objective  function  is  represented  as 


Fk  =  mm(42ax  +  4Lx) 

with  the  constraint  >  0. 


(65) 
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(ii)  The  variable  of  the  object  function  is  A0k-  Then,  considering  the  angle 


+  Aft 


(66) 


and  using  the  equation  of  meshing  with  /3fc,  we  can  determine  the  new  characteristic,  the  piece 

of  envelope  and  the  new  deviations.  Iterations  are  required  to  provide  the  sought-for 

objective  function.  The  final  correction  of  angle  we  designate  as 

Note  1;  The  new  contact  line  (determined  with  4^^)  slightly  differs  from  the  real 

contact  line  since  the  derivative  but  not  is  used  for  determination  .  However, 

as 


close  to  the  real  contact  line. 

Step  3:  The  procedure  discussed  must  be  performed  for  the  set  of  pieces  of  surfaces  Eg  with 
the  characteristic  Lgk  for  each  surface  piece.  Remember  that  the  deviations  for  the  whole 
surface  must  satisfy  the  inequality  >  0.  The  procedure  of  optimization  is  illustrated  with 
the  flowchart  shown  in  figure  21. 


Curvatures  of  Ground  Surface  Eg 

The  direct  determination  of  curvatures  of  Eg  by  using  surface  Eg  equations  is  a  com¬ 
plicated  problem.  The  solution  to  this  problem  can  be  substantially  simplified  using  the 
following  approach  proposed  by  the  authors:  (i)  the  normal  curvatures  and  surface  torsions 
(geodesic  torsions)  of  surfaces  Ep  and  Eg  are  equal  along  line  L-m,  respectively,  (ii)  the  nor¬ 
mal  curvatures  and  surface  torsions  of  surfaces  E^  and  Eg  are  equal  along  line  Lg.  This 
permits  the  derivation  of  four  equations  that  represent  the  principal  curvatures  of  surface 
Eg  in  terms  of  normal  curvatures  and  surface  torsions  of  Ep  and  Et.  However,  only  three  of 
these  equations  are  independent  (see  below). 

Further  derivations  are  based  on  the  following  equations: 


ki  cos^  q  -i-  kn  sin^  q  =  +  ku)  +  -^{ki  -  hi)  cos  2q 


(67) 
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t  =  0.5(fc//  —  ki)  sin  2?  (68) 

Here:  ki  and  kjj  are  the  surface  principal  curvatures,  angle  q  is  formed  by  unit  vectors 
6/  and  e  that  is  measured  counterclockwise  from  ej  and  e;  ej  is  the  principal  direction  with 
principal  curvature  kj]  e  is  the  unit  vector  for  the  direction  where  the  normal  curvature  is 
considered;  t  is  the  surface  torsion  for  the  direction  represented  by  e. 

Equation  (67)  is  known  as  the  Euler  equation.  Equation  (68)  is  known  in  the  differential 
geometry  as  the  Bonnet- German  equation  (see  section  5). 

The  determination  of  the  principal  curvatures  and  principal  directions  for  is  based  on 
the  following  computational  procedure  (see  section  5): 

Step  1:  Determination  of  k^'>  and  for  surface  at  the  direction  determined  by  the 
tangent  to  Lm- 

The  determination  is  based  on  equations  (67)  and  (68)  applied  for  surface  Sp.  Recall 
that  Sp  and  have  the  same  values  of  k^^^  and  along  the  above  mentioned  direction. 
Step  2:  Determination  of  kj^^^  and 

The  designations  kj^^^  and  indicate  the  normal  curvatures  of  S^  and  the  surface  torsion 
along  the  tangent  to  L^.  Recall  that  kj^^^  and  are  the  same  for  St  and  Sp  along  Lp.  We 
determine  and  for  surface  St  using  equations  (67)  and  (68),  respectively. 

Step  3:  We  consider  at  this  stage  of  computation  that  for  surface  Sp  are  known:  kj^^^  and 
k!;^'*  and  for  two  directions  with  tangents  Ti  and  T2  that  form  the  known  angle  fi 
(fig.  22).  Our  goal  is  to  determine  angle  qi  (or  92)  for  fh®  principal  direction  e^®^  and  the 
principal  curvatures  kj^'^  and  k\^j  (figure  22). 

Using  equations  (67)  and  (68),  we  can  prove  that  k!;^^  and  (z  =  1,2)  given  for  two 
directions  represented  by  Ti  and  T2  are  related  with  the  following  equation 

-1- 

A/Tt  — 
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step  4:  Using  equations  (67)  and  (68),  we  can  derive  the  following  three  equations  for 
determination  of  qi,  and 


tan2$i  -  ^(2)_i(i)cos2/i 

(70) 

=  fcW  _  fW  tan  qi 

(71) 

(72) 

Equation  (70)  provides  two  solutions  for  qi  (9^^^  =  +  90°)  and  both  are  correct.  We 

choose  the  solution  with  the  smaller  value  of  qi. 

Numerical  Example  1:  Grinding  of  Archimedes’  Worm  Surface 

The  worm  surface  shown  in  figure  23  is  a  ruled  undeveloped  surface  formed  by  the  screw 
motion  of  straight  line  A7V  {\KN\  =  Up).  The  screw  motion  is  performed  in  coordinate 
system  Sp  (figure  23(b)).  The  to  be  ground  surface  Sp  is  represented  in  Sp  as 

Tp  =  Up  cos  a  cos  6p  ip  +  Up  cos  a  sin  9p  jp  +  {p9p  -  Up  sin  a)  kp  (73) 

where  Up  and  9p  are  the  surface  parameters. 

The  surface  unit  normal  is 


Np  dvp  dvp 

~  |Np|  ’  “  dup  d9p 

Thus; 


(74) 
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np  = 


(u2+p2)0.S 


p  sin  6p  +  Up  sin  a  cos  6p 
—p  cos  6p  +  Up  sin  a  sin  6p 


Up  cos  a 


(provided  cosa^O)  (75) 


As  an  example  the  following  data  will  be  used: 

Number  of  threads,  Ni,  =2 


1 


Axial  diametral  pitch,  Par,  =  8  (— ) 

in 

a  =  20° 

The  radius  of  the  pitch  cylinder  =  1.125  (in) 

(i)  The  screw  parameter  is 


p  =  ^  =  0.125  in. 


(ii)  The  lead  angle  is 


tan  Ap 


L 


0.125 
1.25  ’ 


Ap  =  5.7106° 


The  mean  line  is  determined  as 


^^p(^m,^p)  ,  UfYi 


Tv  + 


+  kp  - 


1.25 


0.125 


^v- 


2  cos  o; 


cos  a 


=  1.3136  in. 


1  1.25 

where  and  —  determine  the  addendum  and  dedendum  of  the  worm. 

iax  ^ax 

The  worm  is  ground  by  a  cone  with  the  apex  angle  7f  =  30°,  and  outside  diameter  8  in.. 
The  inside  angle  =  —88.0121°  provides  the  coincidence  of  both  generatrices  of  the 
cone  and  the  Archimedes’  worm.  The  maximal  deviation  of  the  ground  surface  Sp  from  the 
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ideal  surface  Sp  with  the  above  value  of  is  3  microns. 

The  optimal  angle  =  -94.6788°  has  been  determined  by  the  optimization  method 
developed.  The  deviations  of  the  ground  surface  Sp  from  Sp  with  the  optimal  are 

positive  and  the  maximal  deviation  has  been  reduced  to  0.35  microns  (figure  24). 

5.  Condensed  Information  About  Surface  Curvature 

The  contents  of  this  section  provide  a  condensed  overview  of  about  the  basic  equations 
of  surface  curvatures.  For  further  explanation  of  the  details  please  refer  to  the  books  by 
Nutbourne  and  Martin[6],  Favard[7]  and  Litvin[2]. 

Osculating  Plane 

Figure  25  shows  spatial  curve  L1ML2.  The  osculating  plane  is  the  limiting  position  of 
such  a  plane  that  passes  through  curve  points  Mi,  M,  and  M2  as  Mi  and  M2  approach  M . 

The  osculating  plane  for  a  curve  at  its  regular  point  M  is  formed  by  the  tangent  to  the 
curve  and  the  acceleration  vector  for  the  same  point. 

The  osculating  plane  and  the  curve  are  in  tangency  of  second  order.  The  osculating  plane 
is  an  exceptional  tangent  plane.  The  deviations  of  the  curve  from  the  osculating  plane  are 
of  different  signs  on  the  two  sides  from  the  point  of  tangency,  and  the  curve  is  above  and 
below  the  plane  (see  points  Li  and  L2  in  figure  25).  An  exception  is  the  case  when  the  point 
of  tangency  is  a  rectification  point  at  which  the  second  derivative  r**  of  a  curve  represented 
by  r(s)  is  equal  to  zero.  Here:  s  is  the  arc  length  of  the  curve. 

Space  Curve  and  Surface  Trihedron 

Henceforth,  we  will  consider  two  trihedrons,  the  space  curve  trihedron  and  the  surface 
trihedron.  Each  of  the  trihedrons  is  right-handed,  formed  by  three  mutually  perpendicular 
vectors.  The  concept  of  space  curve  trihedron  is  discussed  when  a  space  curve  is  considered 
in  the  3D  space  and  the  curve  is  not  related  to  a  surface.  The  concept  of  surface  trihedron 
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and  space  curve  trihedron  are  considered  simultaneously  when  the  space  curve  belongs  to 
a  certain  surface  and  the  curve  inherits  some  of  the  properties  of  the  surface  to  which  it 
belongs. 

Space  Curve  Trihedron 

We  consider  a  coordinate  system  that  is  rigidly  connected  to  the  curve.  Position  vector 
OC  =  r(s)  determines  the  current  point  C  of  the  curve  (figure  25);  5  =MC  is  the  length  of 
the  curve  arc;  M  is  the  starting  point. 

Consider  that  a  small  piece  of  curve  L1ML2  is  located  in  the  osculating  plane  no(figure 
25).  Plane  is  perpendicular  to  plane  Hq  and  passes  through  point  M  of  the  curve. 

We  define  the  normal  N  to  the  curve  as  a  vector  that  is  perpendicular  to  the  tangent 
to  the  curve.  There  is  an  infinite  number  of  normals  N  to  the  curve  at  its  point  M.  All  of 
normals  N  belong  to  plane  11;^/  since  the  unit  tangent  t  is  perpendicular  to  TTat.  For  instance, 
vector  N,-  is  one  of  the  set  of  curve  normals  (figure  25).  Two  normals  of  the  set  of  normals 
must  be  specified: 

(i)  the  principal  normal  with  the  unit  vector  m  that  lies  in  the  osculating  plane  IIo  and  is 
the  line  of  intersection  of  planes  IIo  and  Il^r  (figure  25)  and 

(ii)  the  binormal  b  that  is  perpendicular  to  t  and  m  simultaneously. 

We  may  identify  at  each  current  point  of  the  curve  three  mutually  orthogonal  vectors 
(figure  25):  the  tangent  vector  t,  the  principal  normal  m,  and  the  binormal  b.  The  ori¬ 
entation  of  these  vectors  in  a  fixed  coordinate  system  is  varied,  depending  on  the  location 
of  the  point  on  the  curve.  We  may  consider  now  a  trihedron  Sc  as  a  rigid  body  with  three 
mutually  perpendicular  vectors  ec(ic,  jc  ,kc)  that  form  a  right  trihedron  (figure  26).  The 
origin  of  the  trihedron  moves  along  the  curve,  and  the  unit  vectors  ij,  jc,  kc  represent  t,  m, 
b,  respectively.  Unit  vectors  t,  m,  b  are  taken  at  the  current  point  of  the  curve  where  the 
origin  of  trihedron  Sc  is  located  at  this  instant. 

The  representation  of  unit  vectors  t,  m,  and  b  in  terms  of  derivatives  of  vector  function 
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r(5)  is  based  on  the  following  consideration: 

(i)  Unit  vectors  t,  m,  and  b  form  a  right-hand  trihedron  (figures  25  and  26).  Thus 


t  =  mxb,  m  =  bxt,  b  =  txin 

(ii)  Unit  vector  t  is  directed  along  the  tangent  to  the  curve  and  therefore 


(76) 


,  .  dr 
t(.)  =  -  =  r. 


Vector  r,  is  a  unit  vector  since  |  dr  |=  ds. 

(iii)  The  principal  normal  to  the  curve  is  perpendicular  to  the  curve  tangent  t  =  Fj.  The 

derivative  Fj^  =  —  (r,)  is  perpendicular  to  r,,  hes  in  the  osculating  plane  and  therefore  the 
ds 

unit  vector  m  of  the  principal  normal  is  represented  as 


m(s)  = 


r„ 


(iv)  Taking  into  account  the  expression  for  b  in  equations  (76),  we  obtain  the  following 
equation  for  the  binormal 


b(5)  =  t  X  m  = 


r*  X  r„ 


Frenet-Serret  Equations 

The  motion  of  the  trihedren  along  a  spatial  curve  can  be  represented  in  two  components: 
(i)  as  a  translational  motion  along  the  curve  (the  origin  of  the  trihedron  moves  along  the  curve 
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and  the  unit  vectors  of  the  trihedron  keep  their  original  orientation),  (ii)  and  as  a  rotational 
motion  (the  trihedron  is  rotated  as  a  rigid  body  (to  be  coincided  with  the  principal  normal 
nic  and  the  tangent  tc  to  the  curve  at  the  curve  neighboring  point). 

Consider  that  the  origin  of  curve  trihedron  coincides  with  point  M  of  the  curve  and  the 
unit  vectors  t^,  rtic  and  determine  the  instantaneous  orientation  of  the  trihedron(figure 
26).  The  neighboring  point  of  the  curve  is  N  and  \MN\  =  ds, where  s  is  the  arc  length  of 
the  curve.  The  unit  vectors  of  the  trihedron  at  N  are  determined  as  (t*,  m*  and  b*),  where 

t*  =  t,  +  t,cds,  m*  =  rrie  +  m,cds,  b*  =  b^  +  h,Js  (77) 

Here: 

dtc  dnic  -  dbc 

—  "3  1  —  “3  ,  bjj  —  — 

as  as  as 


that  are  taken  at  point  M. 

Frenet-Serret  equations  define  t^c,  rrijc  and  bjc  as  follows  (see  References  [6],  [7]  and  [2]): 


u 

1 _ 

KolTlc 

o 

o 

O 

tc 

= 

rbc  -  /Cote 

= 

—  /Co  0  r 

ntic 

b,. 

— rnic 

0  — T  0 

be 

(79) 


where  Kq  and  r  are  the  curvature  and  torsion  of  the  space  curve  at  point  M.  It  is  evident 
that  in  the  case  of  a  planar  curve,  the  unit  vector  b,.  is  perpendicular  to  the  plane  where  the 
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curve  is  located,  h,c  is  equal  to  zero  and  the  curve  torsion  r  of  a  planar  curve  is  equal  to 
zero. 

Equations  of  and  r  for  a  Parametric  Spatial  Curve 

Consider  that  the  spatial  curve  is  represented  by  vector  function  r(^).  After  derivations 
we  obtain  (see  References  [2],  [6],  and  [7]) 


Tee  ■  m 


Te  X  re$ 


[i^eVee  —  +  {^ezee  —  ^eezeY  +  {yezge  —  VeezeYY^^ 

{xl  +  yj  +  zlf^^ 


The  curvature  obtained  from  equation  (80)  is  always  positive  because  the  principal 
normal  m^  is  located  in  the  osculating  plane  and  is  directed  to  the  center  of  curve  curvature. 
The  curve  curvature  can  be  also  represented  in  the  form 


(81) 


Here:  and  a,,  are  the  velocity  and  acceleration  of  a  point  in  its  motion  along  the  curve 

and  are  represented  as  follows 

V,  =  r/1  (82) 

a.  =  r.4)%r.(^)  (83) 


Obviously,  the  curvature  Ko  can  be  also  represented  as 
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Tee  ■  m 


(84) 


Ko  = 


r 


2 

e 


The  space  curve  torsion  r  is  represented  by  the  equation 

_  (rg  X  rgg)  •  rgflg 
^  "  (r,  X  ree)^ 

In  the  case  of  a  planar  curve,  we  have  (r^  x  rgg)  •  Tgse  =  0  and  r  =  0. 

Surface  Curve  Trihedron 

Consider  a  regular  surface  S  that  is  represented  by 

t(u,  d)  ,  r„  X  re  ^  0  ,  {u,B)  e  A 


(85) 


(86) 


A  curve  on  E  is  determined  if  in  vector  function  r(u,  6)  surface  parameters  are  related 
with  the  equation 


(87) 


Figure  27  shows  two  curves,  Ln  and  that  pass  through  the  same  surface  point  M  and 
have  the  same  tangent.  Curve  is  a  planar  curve  obtained  by  intersection  of  the  surface 
by  the  surface  normal  plane  that  is  drawn  through  the  unit  tangent  t  and  the  surface  unit 
normal  n.  Curve  Lo  is  a  spatial  curve  identified  locally  with  the  orientation  of  osculating 
plane,  the  curvature  and  the  torsion  of  the  curve.  Considering  that  a  spatial  curve  belongs 
to  a  surface,  we  may  determine  more  parameters  for  the  local  identification  of  the  curve. 
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We  have  introduced  in  above  the  curve  trihedron  iS'e(ic,  jo  kg)  where  ic  —  t  is  the  unit 
tangent,  =  m  is  the  curve  principal  normal,  and  kc  =  b  is  the  curve  binormal  (figures 
26  and  27(b)).  In  addition,  we  set  up  now  the  surface  trihedron  5/(1/,  j/,  k/)  shown  in 
figure  27(b).  Here:  i/  =  t  is  the  unit  tangent  to  the  spatial  curve,  j/  =  d  is  the  unit  vector 
that  is  perpendicular  to  t  and  lies  in  the  plane  tangent  to  the  surface  at  point  Af ,  k^  n  is 
the  surface  unit  normal.  Subscript  “/”  indicates  that  the  surface  trihedron  and  its  axes  are 
considered. 

The  unit  tangent  iy  =  =  t  is  determined  as 

t  =  -^,  T  =  r^  +  ref-  =  T^-re^  (fe^O)  (88) 

I  T  I  du  fe 

The  surface  unit  normal  is  represented  as 

n  =  ,  N  =  X  re  (n  =  k/)  (89) 

Changing  the  order  in  the  cross  product  in  equation  (89),  we  can  change  the  direction 
of  n  for  the  opposite  one,  and  provide  5  <  90°,  where  6  is  formed  by  n  and  m.  We  remind 
that  the  direction  of  m  is  the  same  as  (assuming  that  the  curve  is  represented  by  r(s)) 
and  cannot  be  chosen  arbitrarily.  Unit  vectors  t,  d,  and  n  form  the  right  trihedron  Sf,  the 
surface  trihedron. 

Bonnet-Kovalevski  Equations 

Figure.  27(b)  shows  the  curve  and  surface  trihedrons  whose  common  origin  is  located 
at  the  current  point  M  of  spatial  curve  L<,.  Consider  now  that  the  common  origin  of  both 
trihedrons  is  moved  along  Lo  to  the  neighboring  point  N.  Both  trihedrons  will  keep  the 
tangent  t*  to  Lo  at  point  N  as  their  common  axis,  but  one  of  the  trihedrons  will  be  turned 
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with  respect  to  the  other  one  since  the  motion  along  will  be  accompanied  with  the  change 
of  angle  S  formed  by  vector  m  and  n.  Obviously,  the  unit  vectors  of  the  surface  trihedron 
will  change  at  N  their  orientation  with  respect  to  the  orientation  at  M.  Designating  the  unit 
vectors  at  N  by  t*,  d*  and  n*,  we  have 


t"  =  t{s)  +  tsds,  d*  =  d(s)  +  d,ds,  n’ =  n(s)  +  (90) 


where 


(91) 


Bonnet-Kovalevski  equations  express  the  derivatives  t,,  d,  and  in  terms  of  Kg,  and 
t  as  follows  (see  References  [6]  and  [2]). 


t,  =  K^d  +  /C„n  =  Kgjf  +  Knkf 


ds  = -Kgt  +  tn  = —Kgif  +  tkf  y 


(92) 


n5  =  -Knt  -td-  -K„i/  -  tjf 


Here:  Kg  and  t  are  the  surface  normal  curvature,  geodesic  curvature,  and  the  surface 

torsion,  respectively.  The  concept  of  surface  normal  and  principal  curvatures  is  discussed  in 
many  books  on  differential  geometry,  but  the  determination  and  concept  of  Kg  and  t  requires 
additional  explanation  that  is  presented  next  in  this  report. 

Geodesic  Curvature 

Frenet-Serret  equations  (92)  yield  that 
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tj  =  /Com 


(93) 


where 


/c  is  the  curvature  of  a  spatial  curve;  the  curvature  center  lies  in  the  osculating 


plane.  Equations  (92)  and  (93)  yield  that 


/Com  =  Kgd  +  /c„n 


(94) 


Equation  (94)  can  be  interpreted  follows: 

(1)  Figure  29  shows  a  spatial  curve  ou  surface  S.  Unit  vectors  t,  d,  aud  n  represent  the 
surface  trihedron  (figures  29  and  27(b)).  Here:  t  is  the  tangent  to  curve  i,;  d  lies  in  the 
tangent  plane  and  is  perpendicular  to  t;  n  is  the  surface  unit  nornaal.  Unit  vector  m  is  the 

principal  normal  to  Lo  and  lies  in  the  osculating  plane.  Vector  r„  =  K,m. 

(2)  Consider  now  that  the  spatial  curve  L,  is  projected  on  the  tangent  plane  T  and  normal 
plane  N,  respectively.  The  projections  are  designated  by  Lr  and  ijv.  We  emphasize  that 
there  is  no  difference  between  iw  (figure  29(b))  and  i„  (figure  28)  it  they  are  considered 
locally.  Both  curves  have  the  same  normal  curvature  at  the  point  of  tangency  M. 

(3)  Vector  /c„m  is  represented  as  the  sum  of  two  vectors:  a,d  and  «:,n.  The  scalar  a, 
represents  the  curvature  of  curve  Lr,  and  the  scalar  a„  represents  the  curvature  of  curve  i„. 

(4)  Equation  (94)  yields  two  relations 

/Co(m  •  n)  =  /Co  cos  ^  =  /c„  (9^) 

•  d  = /CoSin(5  (96) 
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where  6  is  the  angle  formed  by  vectors  m  and  n  that  determines  the  orientation  of  osculating 
plane  with  respect  to  the  normal  plane.  Equations  (95)  and  (96)  relate  curvatures  Kq  and 
and  angle  8. 

The  direct  determination  of  geodesic  curvature  of  a  spatial  curve  represented  on  a  surface 
is  based  on  the  following  equations: 

(i)  Consider  that  the  surface  is  represented  by  the  vector  function  r(u,^). 

(ii)  A  spatial  curve  is  represented  on  the  surface  as 

T(u(e),e)  (97) 


where  u{6)  is  the  known  function. 

(iii)  The  tangent  to  the  curve  is  represented  as 


^  du 
T  =  r„-  +  r, 


(98) 


(iv)  The  unit  normal  to  the  surface  is  represented  as 


n  =  ,  N  =  r„  X  re 


(99) 


(v)  An  auxiliary  parameter  a  is  represented  as 


.du.o  ^  du 
a  =  r„:.(— )  +2r^9—  +  Tee 


(100) 


(vi)  The  final  expression  for  Kg  is 
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(101) 


K 


3 


T  •  (a  X  n)—  1  N 

TtF" 


(Pu 

dp 


Surface  Torsion 

The  surface  torsion  t  can  be  represented  by  the  equation 

l  =  T  +  «,  =  T+^(«)  (102) 

Thus,  the  surface  torsion  depends  on  the  torsion  r  of  the  spatial  curve  along  which  the 
origin  of  two  trihedrons  is  moved,  and  on  the  derivative  6^  where  8  is  the  angle  formed  by 
the  unit  vectors  m  and  n  of  the  trihedrons. 

The  geometric  interpretation  of  the  surface  torsion  may  be  based  on  the  concept  of  the 
geodesic  line  (see  References  [6]  and  [2]).  A  spatial  line  on  the  surface  is  the  geodesic  one  if 
the  principal  normal  m  at  any  curve  point  M  coincides  with  the  surface  normal  at  M.  The 
geodesic  curvature  of  the  geodesic  line  at  any  curve  point  is  equal  to  zero. 

It  was  proven  in  differential  geometry  that  the  surface  torsion  t  is  the  curve  torsion  of 

the  geodesic  line. 

A  simple  method  for  computation  of  the  surface  torsion  is  based  on  the  equation  that 
has  been  proposed  by  Sophia  Germain  and  Bonnet  (see  References  [6],  [2]  and  [7]).  This 
equation  is 


t  =  0.5(k//  —  Kj)  sin2g 


(103) 


Here:  kj  and  /cj/  are  the  principal  curvatures  of  the  surface  at  point  M  on  the  principal 
directions  with  the  unit  vectors  ej  and  ejj(fig.  30);  q  is  the  angle  formed  by  e/  and  t. 
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Using  equation  (103)  and  Euler’s  equations  that  relate  the  principal  curvatures  and  nor¬ 
mal  curvatures,  we  may  determine  the  solutions  to  the  following  two  problems  (Ref.  [2]); 
Problem  1:  Consider  that  two  directions  in  the  tangent  plane  determined  with  unit  vectors 
and  are  given  (figure  31).  The  angle  n  formed  by  and  is  known.  Also  the 
following  are  known:  (i)  the  normal  curvatures  and  for  directions  of  and 
and  (ii)  the  surface  torsion  given  for  direction 

The  goal  is  to  determine  the  principal  curvatures  acj  and  zc//  for  directions  of  and 
and  angle  (  or  g(^)). 

The  solution  to  this  problem  is  represented  by  the  following  equations  [2] 


(104) 

Kj  =  tan 

(105) 

Kjj  =  -f  cot 

(106) 

Problem  2:  Consider  as  given  t^^^  t^^^  (figure  31),  k^'^  and 

surface  torsions  for  directions  of  t^^^  and  t^^^. 

The  goal  is  to  relate  the 

The  sought-for  relation  is  represented  by  the  equation  [2] 

^(2)  (1)  cot  /X  _  0 

fvn  fvn 

(107) 

where  ji  is  the  angle  formed  by  and 

Numerical  Example  2:  Determination  of  the  geodesic  curvature  Kg  and  surface 
torsion  t  of  Archimede’s  worm  surface 
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Archimede’s  worm  surface  is  represented  as 


r  = 


u  cos  or 


cos  9 


u  cos  Of  sin  0 


pB  —  u  sin  oc 


(108) 


Here  u  and  9  are  the  Gaussian  coordinates,  a  is  the  pressure  angle  and  p  is  the  screw 
parameter.  A  helix  on  the  Archimede’s  worm  surface  is  a  spatial  curve  obtained  by  inter¬ 
section  of  the  worm  surface  by  a  cylinder  of  radius  r,'.  Our  goal  is  to  determine  the  geodesic 
curvature  Kg  of  the  helix  and  the  surface  torsion  of  the  Archemede’s  worm  surface. 

(1)  Geodesic  curvature  Kg 

Taking  into  account  the  -|-  rj  =  r?,  we  can  represent  the  helix  as  follows 


r,-  cos  9 


ri 

u  = - , 

cos  a. 


r(0)  = 


r,-  sin  9 


p9  —  ri  tan  a 


(109) 


The  tangent  to  the  helix  can  be  represented  as 


T  =  ^  =  {-ri  sin  9  ri  cos  0  pf 


(110) 


The  unit  surface  normal  is  represented  by  equations  (99)  and  (75) 


47 


p  sin  0  +  u  sin  a  cos  6 


n  = 


+  V? 


—p  COS  ^  +  u  sin  a  sin  9 


u  COS  a 


(111) 


The  auxiliary  vector  a  is  (see  equation  (100)) 


a  =  (— r,cos0  —  r,sin0  0)^ 


(112) 


Equations  (101)  and  (110)  to  (112)  yield  the  following  expression  for  the  geodesic  curva¬ 
ture 


K 


9 


T  •  (a  X  n) 

|TP 


r,  cos  a 

\/ (P^  +  ’'?)(p^  cos^  a +  rf) 


(113) 


(2)  Surface  torsion  t 

From  Reference  [2]  (F.  L.  Litvin,  1993),  the  principal  curvatures  and  principal  directions 
at  a  surface  point  can  be  represented  by  the  following  equations: 


Ki  = 


e,  = 


Lhi  +  M 
Ehi  -f  F  ’ 

(114) 

+  ve 
ku/ii  +  rel 

{i  =  /,  II) 

(115) 

The  coefficients  and  the  partial  derivative  in  the  case  of  the  Archimedes’  worm  can  be 
expressed  as  follows: 
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L  =  0,  M  = 


pcosa 


(u^  +p2)2 


F  = —psina,  E  =  l 


(116) 


—u^  sin.  a  ±  {u^  sin^  a  +  Ap^u^  +  Ap^)  2  ({=.  I  JJ) 

’■=  2p  ,  ,  ) 


(117) 


Fu  =  (cos  a  cos  6  cos  o;  sin  0  —  sin  a)^ 


(118) 


r/i - —  ^  ucosacos^  p)^ 

\Jp^  +  t£^  cos^  a'' 


(119) 


Angle  q  that  is  formed  by  tangent  T  and  ej  is 


9  = 


cos“^( 


e/-T 

|T| 


) 


(120) 


Considering  that  /c/,  Kn  and  angle  q  are  given,  we  can  obtain  the  surface  torsion  along 


the  tangent  T  as 


t  =  0.5(kj/  -  Ki)  sin  2q 


(121) 


(3)  Computation  results 

The  to  be  computed  point  is  located  on  the  helix  that  belongs  to  the  pitch  cylinder  of 
the  worm.  The  z-coordinate  of  the  helix  point  is  equal  to  zero,  and  the  Gaussian  coordinates 

are 


=  ■,  0  =  — tana 

cos  a  p 
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The  design  parameters  are  the  same  as  in  Numerical  Example  1,  i.e. 


Tp  =  1.25m,  p  =  0.125zn,  o;  =  20° 

The  results  of  computation  are 

/c,  =  0.8257^,  t  = -0.1540— ,  /c/ = -0.3283— ,  k// =  0.0227^ 

in  in  in  in 

6.  Conclusion 

From  the  analytical  study  presented  in  this  report  the  following  conclusions  can  be  drawn: 

(1)  The  kinematics  of  two  CNC  machines  with  6  degrees-of-freedom  has  been  described. 

(2)  The  preliminary  results  of  investigation  by  TCA  of  the  sensitivity  of  helical  gears  and 
worm-gear  drives  to  misalignment  are  represented. 

(3)  A  new  method  for  grinding  of  a  gear  tooth  surface  with  optimal  approximation  to  the 
given  surface  is  proposed. 

(4)  An  algorithm  for  the  execution  of  motions  of  a  CNC  machine  for  the  surface  generation 
has  been  developed. 
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Fig.  2  Transmission  errors  of  helical  gears:  error 
of  crossing  angle  is  A7  =  S.Oarc  —  min 
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Fig.  3  Edge  contact  of  helical  gears:  error  of  pinion  lead 
angle  on  pitch  cylinder  is  AApi  =  S.Oarc  —  min 
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Fig.  4  Transmission  error  of  helical  gears:  error  of  pinion 
lead  angle  on  pitch  cyhnder  is  AApi  =  S.Oarc  —  min 
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Fig.  7  Transmission  errors  of  misaligned  worm-gear  drive 


Fig.  9  Schematic  of  “Phoenix”  machine 
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Fig.  10  Coordinate  systems  applied  to  “Phoenix”  machine 
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Fig.  11  ‘‘Star”  CNC  machine 
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Fig.  12  Pinion  head-cutter 
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Fig.  15  Installment  and  orientation  of  tool  surface 
with  respect  to  ideal  surface  Sp 
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Fig.  16  Tool  surface 
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Fig.  17  Orientation  of  trihedron  St  with  respect  to  5/ 
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Fig.  21  Flow  chart  for  optimization 
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.  22  To  determination  of  principal  directions 
of  generated  surface  Ej, 


Fig.  23  Surface  of  Archimedes’  worm 
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Fig.  24  Deviations  of  the  ground  surface  Sj  from 
ideal  surface  Sp  of  Archimedes’  worm 
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Fig.  25  Osculating  and  normal  planes 
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Fig.  26  Curve  trihedron 
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Fig.  27  Surface  trihedron 
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Fig.  28  Surface  normal  section  and  surface  spatial  curves 
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